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^ , Abstract 

' We study the Next-to-Minimal Supersymmetric Standard Model (NMSSM) as the sim- 

■ plest candidate solution to the ;U-problem in the context of the gauge mediation of super- 

I symmetry breaking (GMSB). We first review various proposals to solve the /U-problem in 

models with the GMSB. We find none of them entirely satisfactory and point out that many 
' of the scenarios still lack quantitative studies, and motivate the NMSSM as the simplest 

, possible solution. We then study the situation in the Minimal Supersymmetric Standard 

I Model (MSSM) with the GMSB and find that an order 10% cancellation is necessary be- 

tween the /i-parameter and the soft SUSY-breaking parameters to correctly reproduce Mz- 
P-i. Unfortunately, the NMSSM does not to give a phenomenologically viable solution to the 

(— I I /i-problem. We present quantitative arguments which apply both for the low-energy and 

high-energy GMSB and prove that the NMSSM does not work for either case. Possible 
. modifications to the NMSSM are then discussed. The NMSSM with additional vector-like 

I quarks works phenomenologically, but requires an order a few percent cancellation among 

' parameters. We point out that this cancellation has the same origin as the cancellation 

required in the MSSM. 
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1 Introduction 



The primary motivation for supersymmetry (SUSY) is to stabilize the smallness of the elec- 
troweak scale against radiative corrections 0, which can be as large as the Planck scale 
if the Higgs bosons are truly elementary. Once the electroweak scale is set in the tree-level La- 
grangian, it only receives logarithmic radiative corrections, and hence its order of magnitude is 
not changed. Moreover, the electroweak symmetry remains unbroken in the Minimal Supersym- 
metric Standard Model (MSSM) in the absence of explicit SUSY-breaking parameters. Therefore, 
one can view the electroweak symmetry breaking as being triggered by the soft SUSY breaking. 
Indeed, the soft SUSY-breaking mass-squared of the Higgs boson can be driven negative due to 
the top quark loop 0] while all the other scalar bosons still have positive mass-squared. In this 
sense, there is nothing special about the Higgs boson. It is just one of many scalar bosons, which 
happens to acquire a negative mass-squared due to the top quark loop. This idea eliminates 
one of the least appealing features of the Standard Model. However, there are at least two open 
questions. First, SUSY by itself does not explain why the electroweak scale is small to begin 
with. Therefore, SUSY makes the smallness of the electroweak scale "technically natural," but 
not truly natural. Second, the MSSM contains one dimensionful parameter (the /x-parameter) , 
allowed by SUSY, in the superpotetial. The natural values of fi are either the Planck mass 
(the only natural dimensionful parameter available) or zero, but recent experimental constraints 
imposed by LEP2 imply that a nonzero /i ^ 50 GeV is required 0. 

SUSY, fortunately, can potentially explain the smallness of the electroweak scale if it is 
broken dynamically 0. The perturbative non-renormalization theorem forbids the generation 
of a mass scale in the superpotential if it is absent at the tree-level. However, non-perturbative 
effects can violate the non-renormalization theorem, and a mass scale can be generated by a 
dimensional transmutation: Ag^^ ~ MpianckC"^'^ l^"', if an asymptotically free gauge theory 
is responsible for SUSY breaking. There has been major progress in building models of dynamical 
SUSY breaking |p, |^, ||, ^ |T0|, [11], |12|, which became possible with the detailed understanding of 



the non-perturbative dynamics of SUSY gauge theories [T^. Furthermore, the so-called gauge 



mediation of SUSY breaking (GMSB) 0, can generate soft SUSY-breaking parameters in 



the SUSY Standard Model in a phenomenologically desired form. Therefore, there is hope of 
understanding the smallness of the electroweak scale in a truly natural manner. 

However, the other question remains largely unanswered: how can the dimensionful param- 
eters in the superpotential naturally be of the order of the SUSY-breaking parameters? There 
have been extensive discussions on this subject in the literature which we briefly summarize in 
Section 2. Unfortunately, many of the proposed mechanisms rely on either small parameters, ac- 
cidental cancellations, or the absence of interactions allowed by symmetries. We find the current 
situation to be rather unsatisfactory. 

A natural direction to follow is to start with a superpotential which does not contain a di- 
mensionful parameter and hope that the electroweak scale is generated solely due to the soft 
SUSY-breaking parameters. The simplest model which can potentially work along this line is 
the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [|l5l, which replaces the /i- 
parameter by the vacuum expectation value of an electroweak singlet superfield. We revisit 
this possibility with detailed quantitative studies in this paper. Unfortunately, our conclusion 
is negative. The NMSSM by itself does not produce a phenomenologically viable electroweak 
symmetry breaking even if we vary the messenger scale. The major experimental constraints 
include Higgs boson and slepton searches. Certain simple modifications can evade phenomeno- 
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logical constraints, but require a cancellation among parameters accurate to a few percent. We 
present all of these points quantitatively in this paper, and hope that our results prompt further 
investigations in understanding the origin of the //-parameter in models with the GMSB. 

The paper is organized as follows. In the next section, we review the situation of the /i- 
problem in models with the GMSB, and discuss various proposals to explain the origin of the 
/i-parameter. We, however, find none of them entirely satisfactory. Even if we accept one of the 
proposed models, it is still necessary to check whether the generated //-parameter is phenomeno- 
logically allowed. We address this question in Section 3, and find that the currently available 
experimental lower bounds on superparticle masses already require a cancellation of order 10% 
between the /i-parameter and soft SUSY-breaking parameters to reproduce the observed Mz- 
Then, as the major part of our study, we present the quantitative results of the electroweak 
symmetry breaking in the NMSSM with the GMSB in Section 4 and find that there is no phe- 
nomenologically viable parameter set even if we vary the messenger scale from 10^ to 10^^ GeV. 
We study various simple modifications of the NMSSM in Section 5, and find that they either do 
not break electroweak symmetry in a phenomenologically viable manner or require a cancellation 
among parameters of order 1%. We finally conclude in Section 6. 



2 The /i-problem in the GMSB 

In this section, we review the /i-problem in the supersymmetric Standard Model in general, and 
also various attempts to solve it in the context of the GMSB. 

The parameter /i is the only dimensionful quantity present in the superpotential of the MSSM 

W = ^iH^H, + X\^he,H, + Xt^QdjH^ + X■JQ^u,Hu. (2.1) 

Here, Qi, Li, Ui, di, Cj are the matter chiral superfields with the obvious notation, and Hu, 
the Higgs doublets. Note that /t is part of the supersymmetric Lagrangian, and hence its origin 
is, naively, unrelated to the origin of the soft SUSY-breaking terms 

Vsoft = m%jHd\'^ + mjjjHul'^ 

^ra^Q\Qj + rfi^^L\Lj + m^^uluj + ra'-^ d\dj + mf^elej 

-mlH^H, + A'lQd^Ha + + ^fL^e^H^. (2.2) 

Phenomenology, on the other hand, dictates that the values of both /i and the soft SUSY- 
breaking masses should be around the weak scale (100 GeV), if SUSY is to be responsible for 
stabilizing the Higgs mass. Therefore, the important question is how the mechanism of SUSY 
breaking can induce a /i-parameter naturally, at the same order of magnitude as the other soft 
SUSY-breaking parameters in the Lagrangian. 

One popular scenario of SUSY breaking is the so-called "hidden sector" SUSY breaking in 
supergravity (SUGRA) |l^. In hidden sector models, SUSY is broken in the hidden sector by 
some mechanism, such as the Polonyi model [l7], gaugino condensation [|1^, or the O'Rafeartaigh 



model []T9l, and the effects of SUSY breaking are mediated to the fields in the supersymmet- 
ric Standard Model only by interactions suppressed by the Planck scale. It therefore requires 
SUSY breaking at a scale A ~ 10^° GeV if the soft SUSY-breaking masses are generated as 
A^/Mpianck- This class of models is able to generate the appropriate soft SUSY-breaking masses 
and /i-parameter given that the /i-term is forbidden in the supersymmetric limit by appropriate 



2 



symmetries, and arises due to SUSY breaking (see, for example, the Giudice-Masiero mecha- 
nism pO[). Hidden sector models have, on the other hand, to face serious bounds imposed by 



flavor-changing neutral currents (FCNC) W^ - Low-energy constraints such as the smallness of 
K^-K^ mixing require the matrices fn^Q ■, ''^ij"' to have eigenvalues degenerate to a few percent, 
or their eigenvectors to be strongly "aligned" with the eigenvectors of the Yukawa matrices Af^ 
(the same is true for A!"]^). Within the SUGRA framework alone, there is no natural mechanism 
to guarantee the degeneracy or the alignment [^. In this case, flavor symmetries are probably 



necessary to ensure either degeneracy p3[ or alignment |23] and suppress FCNC, and some of the 



models presented are also capable of generating the /i-term through flavor symmetry breaking 



2^ , p5[1 . There is also the possibility that string theory generates degenerate squark masses if, 
for instance, the dilaton fleld provides the dominant contribution to the soft SUSY-breaking 
masses [p6| . 

The gauge mediation of supersymmetry breaking is an alternative mechanism which can 
naturally ensure the degeneracy of squarks masses and therefore suppress the dangerous FCNC 
effects. SUSY is somehow broken (hopefully dynamically via dimensional transmutation to gen- 
erate a large hierarchy), and SUSY-breaking effects are mediated to the flelds in the supersym- 
metric Standard Model by the Standard Model gauge interactions. Mediating SUSY breaking 
via gauge interactions is not a novel idea [0, [1^. It allows for SUSY breaking at a lower scale 
(when compared to SUGRA inspired models) and, because all SUSY-breaking effects are trans- 
mitted by flavor blind interactions (the Standard Model gauge interactions), squarks of different 
families have the same mass. This scheme has attracted a lot of interest after the pioneering 
works by the authors of references 0, 0, §], which showed that one can successfully mediate the 
SUSY-breaking effects via gauge interactions with the help of a so-called "messenger sector." 
Their scheme can easily incorporate dynamical SUSY breaking and can explain the origin of the 
large hierarchy between the Planck (string, grand unifled (GUT)) scale and the weak scale. 

The GMSB itself, however, has nothing to say about the /i-parameter unless one introduces 
extra flelds which couple to the particle content of the MSSM. The //-problem in the GMSB is 
the primary interest of this paper. Many solutions to the /i-problem have been suggested by 
different authors and all of them require the introduction of new flelds and/or interactions. Some 
of these solutions will be reviewed shortly. 

In the original models 0, |^, ||, SUSY is broken dynamically in a so-called SUSY-breaking 
sector and the breaking effects are transmitted to the supersymmetric Standard Model via a 
messenger sector. The energy scale of the messenger sector is given by A ~ lO'^-lO^ GeV. There 
are, however, models which do not have a separate messenger sector so that the sector which 



breaks SUSY dynamically is directly coupled to the Standard Model gauge group 0, |TT], [T2[ . In 
this case, the effective messenger scale tends to be much higher. For our purposes it is enough 
to employ a simple version of the messenger sector, as in the original models, and take the 
messenger scale A free parameter. 

The messenger sector can be described by the superpotential 

Ws = \XS^ + kS<^+<^- + KgSqq + KiSll, (2.3) 

where S" is a singlet superfleld, $^ are charged under a U{1) associated with the SUSY-breaking 
sector and are singlets under the Standard Model S'f/(3) x SU{2) x U{1) gauge group. The 
superfleld q (g) transforms as a (3(3), 1, ±1/3) under the Standard Model, while / (f) transforms 
as (1,2,^1/2). 
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We assume that the scalar components of acquire negative SUSY-breaking masses-squared 
due to its interaction with the SUSY-breaking sector (usually accomplished by the so-called "mes- 
senger U{1)" gauge interaction ||), and the potential associated with the scalar component 
of S reads 

Vs = -|m$p(|$+p + + + IkS^-]"^ + + A^^p, (2.4) 

neglecting terms containing / or q. It is easy to see that the scalar and the F components 
of S acquire vacuum expectation values (VEVs) {S) and (Fs) and therefore q and / acquire 
supersymmetric masses proportional to (S) and SUSY-breaking masses-squared proportional to 
(-^5) -El This effect feeds down to the MSSM through loop corrections. Gauginos acquire Majorana 
masses at one loop, while sfermions acquire SUSY-breaking masses-squared at two loops. The 
calculation of these soft SUSY-breaking parameters was done long ago (see 0, and its result 
is well known. At the messenger scale: 

M, = ^nB, (2.5) 

4:71 

Here and below, all = gf /^'k are in the SU{^) normalization, B = (Fs) / {S) in the messenger 
sector discussed above, and n determines the number of messenger sector superfields responsible 
for mediating SUSY breaking. In the example we described above, which will be referred to as 
the model with the minimal GMSB, n = 1. y is the hypercharge of the particle, C2 = 3/4 for 
weak SU{2) doublets (zero for singlets) and C3 = 4/3 for color triplets (zero for color singlets). 
Eq. ( p.6|) guarantees that squarks of different families are degenerate at the messenger scale and 
therefore FCNC effects are safely suppressed. It is interesting to note that, for small n, gaugino 
masses and sfermion masses are comparable. For very large n, on the other hand, sfermion 
masses can be significantly smaller than gaugino masses (by a factor ^/n). 

In the mechanism described above, trilinear couplings are not generated at the same order 
(in loop expansion) at the messenger scale. This is not the case in general, and some models can 
generate trilinear couplings with values comparable to the other soft SUSY-breaking parameters 
even at the messenger scale [|TT|. We will, for most of our discussions, consider 



^;f„,.(A) = 0, (2.7) 

unless otherwise noted. 

The GMSB does not generate a /x-term because of the non-renormalization theorem. There- 
fore /i is an input of the model, and, because it has dimensions of mass, its only nonzero natural 
value is Mpianck {Mstring, Mqut) ■ This is clearly not allowed phenomenologically. The /z-term 
must, therefore, be forbidden at the Planck scale (by, say, a Z3 symmetry) and generated dynam- 
ically. Below, we review various attempts to generate the /x-term in the context of the GMSB. 
The following list is not meant to be exhaustive and our descriptions of the various attempts are 
by no means complete. The review below only intends to show that many attempts have been 
made while none of them appears to be entirely satisfactory. 



^There is a run-away direction q — q, I = I in this potential ||27|| . This problem can be avoided by introducing 
more S fields to the messenger sector. Such details are, however, irrelevant for the rest of our discussion. 
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The simplest solution would be to introduce a term in the superpotential 



W D kSH^H^, (2.8) 

where S is the singlet superfield in Eq. (|2.3| ). In such a scenario /i = k{S) and m\ = k{Fs). ml 
is the SUSY-breaking Higgs mixing mass-squared in Eq. ( p.2| ). 

Phenomenology imposes that both /i and ^Jrn^ are of the order of the weak scale, unless 
one is willing to accept a drastic cancellation among parameters to reproduce the observed Mz- 
Therefore, 

ik{S)y ~ k{Fs) ~ (100 GeV)2, (2.9) 



(S) 

and 



k{S) ~ 100 GeV (2.10) 

"^3 = -"^- (2.11) 

This situation is already excluded experimentally. Eq. ( p^.5| ) states that the gluino mass is 
given by (03/47?) (F^)/ {S), and if Eq. ( p.lO| ) is satisfied one would arrive at Mg ~ 1 GeV, which 
is unacceptable. The same is true for all the other soft SUSY-breaking masses. This is a general 
consequence of Eq. ( |2.11|) . It implies that ^Jm^ ^ /i if all experimental bounds on the SUSY 
spectrum are to be satisfied, while SU{2) x f/(l) breaking requires Eq. ( p.9|) . Some authors refer 
to this puzzle as the /i-problem in the GMSB [ p8[] . 

Another simple solution that does not require the introduction of any extra superfields into 
the theory couples the Higgs superfields to the q superfields present in Eq. ( |2.3| ) P5| . In the 
minimal messenger sector 0, |^, one may have, instead of q and I, a complete 5-1-5 multiplet of 
SU{5) to preserve the gauge coupling unification. One can also use a 10 -|- 10 for this purpose, 
and generate gaugino masses and scalar masses-squared with n = 3. In this case, one can couple 
the components Q in 10 that have the same quantum numbers as left-handed quark doublets 
and the components u that have the same quantum numbers as right-handed up quarks (or their 
corresponding components in 10) to the Higgs doublets. Explicitly, W D XiHdQu + \2HuQu. 
This will induce, in the Lagrangian, a one-loop term proportional to 

16^7^^ s^s ■ ^^-^^^ 

The F vacuum expectation value of S will generate /i ~ je^'^^ ^^"^ "^i — le^ ('^^) ' ^§^i^ 
one runs into Eq. (|2.11|) and must hunt for other solutions. 



All of the models described above couple the MSSM Higgs superfields to those in the mes- 
senger sector. Not only did we encounter the problem of Eq. ( p.llD , but some of the coupling 
constants introduced had to be made fairly small because of the magnitude of {S) and {Fs). 
Another class of solutions tries to get around this issue by introducing another singlet superfield, 
whose vacuum expectation value would generate the /x-term. 

One motivation for such models is to utilize the extra singlet to solve the doublet-triplet Higgs 
splitting in SU{b) grand unified theories via a sliding singlet mechanism [Q. This mechanism 
is known to be unstable against radiative corrections if the soft SUSY-breaking parameters are 
generated at a scale higher than the GUT scale, but can be stable for the low-energy GMSB 



30| . Ciafaloni and Pomarol claim that such a solution would generate a viable /i-term. We 
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believe, however, that the conditions that they impose on the soft SUSY-breaking parameters 
can never be satisfied in the context of the GMSB, where all soft SUSY-breaking masses are 
tightly related. We will comment on this in Section |[ 

The simplest model with the addition of an extra singlet one can imagine, referred to as the 
NMSSM |T^, involves substituting the /i-term in the MSSM superpotential by 



\HaK,N-^N\ (2.13) 

The minimization of the scalar potential for Hd-, Hu and at the weak scale should produce 
VEVs Vd and Vu for both Higgs bosons, thus breaking SU{2) x f/(l), and x for the singlet, /x 
would be equal to Ax. The mg term would arise due to renormalization group (RG) running 
of the A-term XAxEdH^N from the messenger scale to the weak scale. m\ would be equal to 
XAxx. 

Dine and Nelson claim that this model does not work for the low-energy GMSB. A detailed 
analysis was not presented in their paper, and we will explain the problem in Section ^. They 
suggest the introduction of an extra light pair of q' + q' and /' + P as a means to produce a 
viable spectrum. They did not, however, publish a quantitative analysis of the model, and say 



nothing about its naturalness. Agashe and Graesser |^ study this scenario and show that there 
is indeed a solution, but it is fine-tuned. They present a possibility to ease the fine-tuning by 
employing many lepton-like messengers while keeping the number of quark-like messengers small. 
In Sections § and ^, we analyze in great detail the case for both the high- and low-energy GMSB. 

There are ways of giving A^ a VEV which are not related to electroweak symmetry breaking. 
In Ref. two mechanisms are introduced, neither of them very appealing, where the A^ VEV 
is generated at the messenger scale. Namely, 

W D -htsN'^S (2.14) 

or 

W D -kqNqq - kiNll (2.15) 

in addition to the NMSSM. S, q and / are the messenger sector superfields present in Eq. (|2.3| ). 
In the case of Eq. ( |2.14| ), a potential 

Vjv = \kN^ + ks{S)N\^ - ksN^{Fs) (2.16) 

is generated for A^ in the presence of (S) and (Fs) VEVs. If one assumes ks to be small, 

N develops a VEV x = \f^^^^, and fi = ^^Jks{Fs) assuming all other couplings to be of 
order one. It is easy to see, a posteriori, that ks must indeed be small if one is to generate 
a phenomenologically viable fi. Unfortunately this case requires that the soft SUSY-breaking 
masses- squared ~ {ai/4:7i)'^{{Fs) / {S)y and /z^ ~ ks{Fs) are accidentally of the same order of 
magnitude .0 

The superpotential coupling Eq. (|2.15| ) would lead to a potential 



■I- Although ks has to be small, its smallness is natural in the sense of 't Hooft. It can be interpreted as being 
generated due to the breaking of some global symmetry, such as TV ^ e^^^/'^N and -ffi.2 e^'^*/'^i/i_2, while S is 
invariant. This type of symmetry would also explain the suppression of a term N S"^ in the superpotential, which 
would be of order {ksY ■ 
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The linear terms in N arise via tadpole one-loop diagrams involving g's and Vs. This would lead 
to = p-32^^^^|fp — \- (l q)- Again kg and ki would have to be small. This solution still 
faces the problem of explaining why a term NS"^ is not present in the superpotential. Note that 
the presence of such a term would lead to an unacceptably large VEV for A^. One may argue, 
however, that this is "technically natural" because the absence of a term in the superpotential 
is preserved by radiative corrections. An even more serious problem is the need to suppress the 
kinetic mixing / d^OS^N + h.c. to ensure -C Fs; an unacceptably large m| = XF^ would 
be generated otherwise. An order unity kinetic mixing can be induced via radiative corrections 
between the ultraviolet cutoff, say the Planck scale, and the messenger scale, and the bare 
parameter has to be chosen very carefully so that the unwanted mixing term can be canceled 
at the messenger scale. This kinetic mixing can be forbidden if there are two sets of messenger 



fields and if the field couples off-diagonally, e.g., W = Nqiq2 etc |3^. Then the tadpole term 
mentioned above is also forbidden, but a negative mass squared for the N field can be generated 
instead. This would lead to the NMSSM in a successful manner; again the parameters must be 
carefully chosen as in the NMSSM with extra light quark pairs (see Section 5). 

Another solution with extra singlets, which points an interesting way around Eq. (ICTI) , was 



suggested by Dvali, Giudice and Pomarol [g8|. Their idea is to generate the /i-term via the 



following one-loop effective term in the Lagrangian: 



^4^^^^™_^_^^v:^^ (2.18) 

where Da is the supersymmetric covariant derivative. This works because cancels 9'^ in S, 
while leaving 9"^ in S'^ . Then the integral over d9'^ can be done and the /i-term is generated, 
while m| is not. The m\ term would arise at higher loops, or via some other mechanism. 

An explicit realization of this mechanism is the following. Suppose a singlet field N 
acquires a linear term M'^N in the superpotential due to its coupling to the messenger sector. 
Then the superpotential W = N{Y'^ + HuHct — M^) leads to a minimum with (Y) = M, 
N = 0. However, by further coupling to the messenger superfields, i.e. Nqq etc, a one-loop 
diagram of messenger fields generates the operator J d'^9NS'^{S'^S)/{S^S), which contains 
V ~ Y^N{Fs)^/{S). Note that this is the same linear potential generated in the case of 
Eq. (pA5|) . This tadpole term induces a VEV for of order (A^) ~ (Fs)^ / (S) / {Y^ which 
is of the order of the weak scale if (Y)"^ = ~ (-^s)- The Y field plays a crucial role: it slides 
to cancel the F-component VEV of A^ before the tadpole is added and, after SUSY is broken, 
its VEV is shifted and leads to (F^) = m\ ~ /i^, as required by phenomenology. Note that 
the /i-parameter obtained here can be understood as a consequence of the effective Lagrangian 
Eq. ( [2.18| ), which is generated upon integrating out A^ and Y before substituting the effect of 
the VEVs of S. 

The necessary linear term (M'^N) in the superpotential for A^ can be easily generated by 
the kinetic mixing between A^ and S or also by other mechanisms, as pointed out in reference 



34| . One apparent drawback of this realization is that one needs a set of new fields whose 
interactions are arranged in a rather special way. Furthermore one would expect the presence 
of a term proportional to SHuHd in the superpotential. This happens because both 5* and A^ 
couple to the messengers, that is, W D Sqq + Nqq, and have, therefore, the same quantum 
numbers. We have already argued that a couphng SHuHa has to vanish (see Eq. ( p.8|) ). Finally 
we point out that this model also suffers from the cancellation problem present in the MSSM 
(see Section 3). 
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Dine, Nelson, Nir and Shirman H] suggest yet another way of generating a /i-term with the 



introduction of an extra singlet. It was inspired by flavor symmetry models in |2J], and resembles 
a modified version of the NMSSM + Eq. 

W D Xn HdHu + + A„ — (2.19) 

M^; M^; ^ Mil 

where Mpi is the Planck mass. When m = 2, n = 1 and p = 2, it is easy to check that 
~ Xn\f{Fs)- We assume the other couplings to be of order 1. It is also easy to see that one 
would require a very small, carefully chosen coupling A„ in order to guarantee /i ~ 100 GeV. It 
is worth noting that this mechanism does not generate an mg term. 

At last we would like to mention another interesting possibility, pointed out by Yanagida 
and Nilles and Polonsky |]SB|. Their models utilize the accidental equality {Ay)sb/M^^^^^ ~ 



[a/ATr)"^, where Adsb ~ 10^ GeV is the scale of dynamical SUSY breaking (DSB) in models with 
the low-energy GMSB and M^^ = Mpz/a/Stt the reduced Planck mass. By introducing a new 
SUSY-preserving sector with strong gauge dynamics, Yanagida's model generates a VEV for 
the superpotential which cancels the cosmological constant from the DSB sector. The constant 
superpotential in turn generates a //-term of order Adsb(Adsb/^^*)^^^ ~ («/4vr)^ADSB ~ 1 TeV. 
The phenomenology of this model is the same as the previous one (see Eq. (|2.19| )). The model by 
Nilles and Polonsky makes use of the Planck-scale suppressed Kahler potential, / d^9N{z* z) / M^, 
where 2; is a chiral superfield in the DSB sector with an F-component VEV. This operator may 
be present at the tree-level, but may also be generated by gravitational effects. It generates 
a tadpole for the singlet N: V = (A^sb/^*)^- Together with the |iV^ superpotential of the 
NMSSM, it generates a VEV for of order (A^) ~ {A'^^^/M^)f/'^ . Even though these models 
generate the correct /i-term of order the weak scale in the models with the low-energy GMSB, 
this would not work for the high-energy GMSB. 

We consider that none of the mechanisms outlined above are entirely satisfactory. Most 
of them require a very specific choice of parameters and the introduction of extra matter at 
or slightly above the weak scale. Furthermore, most of them have not been studied quantita- 
tively (see, however, Ref. ||3^), and there is no guarantee that they indeed generate the correct 
electroweak symmetry breaking pattern and an experimentally viable spectrum. And last, but 
not least, there is no study of how natural such a solution is, given that a viable pattern of 
electroweak symmetry breaking can be generated. 

It is, therefore, part of our goal to study the simplest of the models mentioned above in detail. 
Before that, we would like to review the status of electroweak symmetry breaking in the MSSM, 
where the /i-term is introduced "by hand." We will point out that, in the case of the GMSB, the 
current lower bounds on superparticle masses already require an order 10% cancellation between 
the /i-parameter and the soft SUSY-breaking parameters. 



3 The /i-parameter in the MSSM 

We reviewed various proposals to generate the /i-parameter in models with the GMSB. In this 
section, we review how electroweak symmetry breaking occurs in the MSSM, assuming that the 
/i-parameter and m\ are somehow generated. In particular, we point out a need for an order 10% 
cancellation between /i-parameter and soft SUSY-breaking parameters in models with the GMSB 
given the current experimental lower bounds on superparticle masses. Note that the case of the 
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NMSSM is different because the /i-parameter is generated together with electroweak symmetry 
brealfing and hence the two problems cannot be clearly separated. This will be discussed in the 
next two sections. 

The tree-level Higgs potential in the MSSM is given by 

V = ml\H,\^ + ml\HX-mliH,H^ + cx.) + ^{HlaHd + HlaH^y + ^i\H^^^ (3.1) 

o o 

where the mass parameters involve both the supersymmetric /x-term and the soft SUSY-breaking 
terms, 

ml = i2^ + m]j^, (3.2) 
ml = ^? + m]j^. (3.3) 

In the MSSM, one can show that the vacuum can always be gauge rotated to the following 
configuration 

The two expectation values need to satisfy f J + f ^ = = (174 GeV)^ in order to reproduce 
the observed M^, and it is conventional to parametrize them by Vd = f cos/3, = f sin/3. The 
minimization condition of the potential can be rewritten in the following form: 

Ml ^ mjr^ - mjr tan^ 6 

-f = + W.fl-i • P'^) 

2ml = (2;.= + ™!,, + mi, ) sill 2/3. (3.6) 

It has been claimed that electroweak symmetry breaking is natural in the MSSM because 
m^^ is easily driven negative due to the presence the top Yukawa coupling in its RG evolution. 
In models with the minimal GMSB such as the original ones in [|^ ^, the boundary condition 
for the supersymmetry breaking parameters are given by Eqs. ( |2.5| , |2.6| , p.7|) . A simple one- 
loop approximation is valid in the case of the low-energy GMSB because of the small logarithm 
between the messenger scale A and the electroweak scale, and one finds 

ml^iMz) ^ ml^{Mz) - ^^^t^^' log (3.7) 

which is always negative. 

The need for a cancellation between the //-parameter and soft SUSY-breaking masses can 
be seen as follows. Experimental constraints bound the superparticle masses from below, which 
hence set a lower limit for the ratio B = {Fs)/{S). Therefore one finds that is bounded 

from below. On the other hand, in order for the observed M| to be reproduced, the /i-parameter 
is constrained by Eq. (|3.5|). For a moderately large tan/? > 2, mjj- can be completely neglected 
and one finds 

2 2 fo o\ 

— fi -ruH^. (3.8) 

This equation requires a cancellation between /t^ and (negative) m^^ to reproduce M\l2 ~ 
(70 GeV)^ correctly. The degree of cancellation is given by (M|/2)//t^.| 

^The degree of cancellation is defined as follows: it is a percentage quantity that measures how much a given 
input parameter (in this case ^^) is free to vary before a given output parameter (in this case M^) changes 
significantly. Explicitly, the degree of cancellation is (d(log Aff )/d(log/i^))^^. This definition corresponds to the 
inverse of the Barbieri-Giudice function ISql. 
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To determine the lower limit on |^«|/^|, we consider a number of experimental constraints 
|39| . One is that the gluino must be heavier than 190 GeV, which becomes stronger if the 
squarks have comparable masses. The second is that the right-handed selectron must be heavier 
than 80 GeV.[] For large tan/?, the right-handed stau may become rather light; we then require 
rrif > 55 GeV if it decays into tau and a neutralino or gravitino, and rrif > 73 GeV if it does 
not decay inside the detector. We also considered the lightest chargino to be heavier than 63 
GeV. The most recent lower bound on the chargino mass |^ is ^ 67 GeV, which leaves our 
analysis virtually unchanged. 

Let us first discuss the case of the minimal low-energy GMSB with small tan (3 to make the 
argument clear. Here we consider A ~ 10^ GeV. The gluino mass constraint requires 

5 > 23 TeV. (3.9) 

This bound itself is independent from the messenger scale. However, the gluino mass bound 
depends on the mass of the squarks, and it strengthens if the squark masses are comparable 
to the gluino mass. For the minimal low-energy GMSB, squarks are significantly heavier than 
the gluino and we can use the bound above. A more stringent constraint is derived from the 
requirement that the right-handed sleptons are heavier than 80 GeV. Including the one-loop 
renormalization group evolution and the D-term, we find 



= 2.89 X 10^^52 -0.232M| COS 2/3. (3.10) 
Therefore we find 

5 > 39 TeV (3.11) 
for the most conservative case cos 2/3 = — 1. With this lower bound we find 

3 y 4:71 J 
Using the one- loop running of mfj^, we obtain 

mjj^{Mz) < -(260 GeVf (3.13) 

and as a result of the minimization condition, 

/i > 250 GeV. (3.14) 

This requires a cancellation of 7% in order to obtain the correct M|. Even though this level 
of cancellation is not of immediate concern, this analysis shows the need for a certain amount 
of cancellation which will become worse as experimental lower bounds on superparticle masses 
improve. 



^This bound depends on the mass of the neutrahno into which the selectron decays. However, since /i turns 
out to be large, it is a posteriori justified to assume that the lightest neutralino is almost pure bino. Then the 
GMSB predicts the relation between selectron and the bino masses, and hence we have a fairly reliable lower 
bound. 
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tan(P) 

Figure 1: Lower bounds on n in models with the GMSB subject to the constraint Mz = 91 GeV 
and to the lower bounds on superparticle masses (see text), (a) as a function of the messenger 
scale, for tan/3 = 2, 10, and 30, and (b) as a function of tan/3 for a fixed messenger scale of 
10* GeV. 



As it is clear from the argument above, the actual lower bound on /i depends on the messenger 
scale and tan/3. We have studied this issue numerically using the experimental bounds quoted 
above and found the lowest possible value of /i as a function of the messenger scale. In Fig. |l|(a) 
we present bounds for three values of tan /3. The lower bound on fi comes from one of the various 
experimental constraints. For instance, the tan j3 = 2 case is dominated by the lower bound on 
mg^ up to a messenger scale of 10^^ GeV, after which the gluino mass bound is more important.^] 
The case of tan/? = 10 has a similar behavior. The situation is more complex and interesting 
for tan/5 = 30. For a messenger scale of up to 10^° GeV, the stau is the lightest supersymmetric 
particle (except for the gravitino). It decays inside the detector to tau and gravitino for the 
lowest messenger scale, but leaves the detector without decaying for higher messenger scales. 
This stable stau provides the strongest constraint. From messenger scales above ~ 10^^ GeV the 
stau decays inside the detector to tau and neutralino. This bound dominates up to ~ 10^^ GeV, 
when the gluino bound dominates. The chargino bound is comparable to that of the gluino for 
the GUT scale {Mgut = 1-86 x lO^*^ GeV). 

In Fig. 0(b) we show the minimum value of /i as a function of tan/? for a fixed messenger 
scale (A = 10^ GeV). The tan/3 dependence can be easily understood as follows. Starting from 

I'We only analyzed the case for one messenger (n=l). For larger n the gluino bound becomes less important 
and the slepton bounds dominate up to the GUT scale. 
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low tan/3, increasing tan/3 decreases the top Yukawa coupling, and hence receives a less 
negative contribution from the top-stop loop. Therefore a lower value of /i is allowed. This 
part is dominated by the cr bound. However beyond tan/3 ~ 20, the bottom and tau Yukawa 
coupling become important. In fact, the scalar tau mass is pushed down both because of the loop 
effect and left-right mixing, and the experimental lower bound on B becomes stronger. Beyond 
tan/3 ~ 30, the stau does not decay inside the detector for this choice of the messenger scale 
and the constraint is even more stringent. This in turn leads to a more negative and hence 
a larger /x. 

Combining both the messenger scale dependence and tan/? dependence, we conclude that 
the most conservative current limit is 



H > 160 GeV. (3.15) 

The required cancellation between /i^ and soft SUSY-breaking parameters in order to reproduce 
the observed Mz is Mf /2/i^ = 16%. Note that this level of cancellation is the absolute minimum, 
and a more accurate cancellation is required for most of the parameter space. 

In the case of minimal supergravity models, where all scalars have the universal SUSY- 
breaking mass-squared ttiq, all gauginos have mass M1/2 and all A-terms are given by AJ = AqXJ 
for f = u,d, I, at the GUT scale, the situation appears to be much better. The renormalization 
group equations can be solved numerically for each choice of tan/3. As an example we take 
tan (3 = 2 and find 

m^^ = ml + 0.50M^/^, (3.16) 
m^^ = -0.32mg - 2.49M2/2 - 0.05^2 - O.2OM1/2A0, (3.17) 
m? = ml + 0.l5M^/2. (3.18) 

By requiring rrig > 80 GeV and M1/2 > 60 GeV (this is a rough bound inferred from the gluino 
bound Mg ^ 190 GeV), we find /x > 82 GeV. This basically does not require any cancellation, 
since M|/2/i2 = 65%. 

The situation can be somewhat ameliorated in the MSSM if there is a Fayet-IUioupoulos 
D-term for the f/(l)y gauge group. Such a D-term is known to arise in many ways, such as 
kinetic mixing of the U{1)y and U{l)mess gauge fields [0. The running of all the parameters 



remains the same except that one adds another contribution from the Dy at the weak scale. If 
the sign is appropriate, it increases m? — rrig + DY and rnfj^ m|^^ -|- |-Dy (less negative) while 
decreasing m|^^ — |-Dy. All of these help push the parameters relevant for electroweak 

symmetry breaking in the right direction. Larger m| reduces the lower bound on (F5) / (S), and a 
less negative m|^^ is also welcome. Therefore the sensitivity to fi (required cancellation between 
H and soft SUSY-breaking parameters) in the MSSM can be improved in the presence of a Dy 
with the appropriate sign. 

We will see in the next two sections that the situation in the NMSSM is much worse. There 
is no phenomenologically viable solution to electroweak symmetry breaking. One can modify 
the model to generate a large negative mass-squared for the singlet field and then find a viable 
solution. This solution also requires a cancellation among parameters which has the same origin 
as the cancellation present in the MSSM. We will also see that the addition of the Fayet- 
lUiopoulos D-term does not improve the situation within the NMSSM. 
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4 The NMSSM with the GMSB 



In this section we study the feasibihty of implementing the GMSB in the framework of the 
Next-to-Minimal Supersymmetric Standard Model (NMSSM). We begin our presentation by 
introducing the NMSSM: its particle content, superpotential, and soft SUSY-breaking terms. 
We briefly review the major steps in our analysis: the boundary conditions for the breaking 
terms, the RG evolution, and the minimization of the weak-scale one-loop effective potential. 
We then describe the results of a numerical scan of a large portion of the model's parameter 
space. We find that it is impossible to evade the present-day experimental constraints. We 
further strengthen this argument by providing a semi-analytical explanation for the inevitability 
of this conclusion. 



4.1 The NMSSM 

The NMSSM represents an attempt to solve the //-problem of the MSSM in the simplest and 
most direct way: the spectrum of the MSSM is augmented by a gauge singlet superfield N, which 
couples to H^Hu and plays the role of the //-term once it develops a nonzero vacuum expectation 
value [0. The original /i-term is banned from the theory so that there are no dimensionful 



parameters left in the superpotential. 

The VEV of the scalar component of is determined by minimizing the scalar potential 
with respect to Hd, Hu, and N simultaneously. It is natural to expect the VEVs to be of the 
same order of magnitude for all three fields, thus generating an effective //-parameter of order 
the weak scale, as required by phenomenology. 

The complex scalar N introduces two additional degrees of freedom to the Higgs sector. 
Therefore, the particle spectrum of the NMSSM contains three CP-even Higgs scalars, two CP- 
odd Higgs scalars, and one charged Higgs scalar. Immediately, there is a problem: one of the 
pseudoscalar Higgs bosons is massless. This happens because the superpotential W = XNH^Hu 
has a Peccei-Quinn symmetry — > Ne'^°' , H^Hu — > HfiHue~^'^. This symmetry is spontaneously 
broken by the VEVs of the fields, making one of the pseudoscalars massless. 

The standard solution to this problem is to introduce a term cubic in A^, which explicitly 
breaks the symmetry mentioned above. This term is allowed by the gauge symmetries of the 
model and does not contain a dimensionful coupling constant, so it is generically expected to be 
present in the superpotential. One, however, still has to worry about a light pseudoscalar Higgs 
boson. As we will show shortly, its mass can also be small because of the presence of a different 
(approximate) U{1) symmetry. 

Overall, the only change made to the MSSM superpotential is the following: 

fiHdH^ ^ XNHdHu - (4.1) 

while the corresponding change to the soft SUSY-breaking part of the potential is: 

- ml{HaH^ + c.c.) -{XA^NH^H^ + hAkN^" + h.c.) + ml\N\\ (4-2) 

One can determine the VEVs of the Higgs fields Hd, Hu, and N by minimizing the scalar 
potential, which at the tree-level consists of the F-terms, Z)-terms, and soft SUSY-breaking 
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terms: 



Vf = \XHdHu - fciVY + A2|iV|2(|iJ,|2 + \Hu\''), 
Vd = ^f{H\aH, + HlaH^r + ^-^{\H,\' - 

Vsoft = mljHdl^ + mljH^l^ + ml\N\^ - {XA^H^H^N + h.c.) - (^^^fciV' + h.c^ . (4.3) 

An important fact to notice is that both Vf and Vd remain unchanged when Hd, Hu, and are 
all rotated by the same phase. In fact, only the soft SUSY-breaking A-terms are not invariant 
under this transformation. This can be potentially dangerous, because we, in general, consider 
the A-terms to be zero at the messenger scale, and their sizes at the weak scale are determined 
by the RG evolution. If the generated values of Ax and A^ are not large enough, our scalar 
potential has an approximate U{1) symmetry. This symmetry is spontaneously broken by the 
vacuum expectation values of the Higgs fields, and, as before, we have to worry about a light 
pseudoscalar Higgs boson. 

We denote the VEVs of the neutral components of the Higgs fields by Vd and Vu, as in 
Section H, and the VEV of the singlet field by x: 

{N) = X. (4.4) 

As a function of these VEVs, the potential has the form 

-{XAxVdVuX + h.c.) - (-AfcX^ + h.c.) + ^ {\vd\^ - \vu\^)^. (4.5) 

o o 

It is well known that some of the Higgs boson masses receive significant contributions from 
radiative corrections. In our numerical analysis we account for that by employing the one-loop 
effective potential 

Vl;:n{vd = V:Zrai{v..^) + ^STrA1^(^,) (log^^^^ - ^) . (4.6) 

In this expression M.'^{vi) is a field-dependent scalar mass-squared matrix, and /i is the MS 
renormalization scale. As we have indicated explicitly, the values of the various parameters 
entering depend on the choice of this scale. To the leading order this dependence is 

canceled when the second term on the right hand side of Eq. ( [4.6| ) is included, and the result 
of minimizing V^~''°°^ is less sensitive to the choice of the scale where one stops running the 
RG equations. (Canceling out this dependence completely would require calculating radiative 
corrections to all orders.) 

The matrix depends on the field VEVs Vi through the Yukawa couplings of the Higgs 
fields to various other particles. What plays a crucial role here is not the absolute values of the 
masses, but rather the rate of their change as one changes Vi. Therefore, the most important 
contribution comes from the field-dependent masses of the top quark and squarks, which have 
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the largest Yukawa coupling. Denoting their mass eigenvalues by rrit, m^^ and m^^ respectively, 
the contribution to V^~^°°p from radiative corrections due to these states is 

AV 



327r2 




-2m^(..)(log^-Jl . (4.7) 



4.2 Numerical Analysis 

In models with the GMSB the values of the soft SUSY-breaking terms are specified at the 
messenger scale by Eqs. ( ^.5] ), ( ^^.6] ) and ( [2.71 ). Their values at the weak scale can be determined 



by solving the RG equations given in [Appendix A . 



The model has five input parameters: ht, A, k, B, and n. (Note that the only dimensionful 
input parameter is B, and its magnitude will determine the overall scale of the VEVs and the soft 
SUSY-breaking masses.) There are, however, two constrains which must be satisfied at the weak 
scale: v = ^ + f ^ = 174 GeV and htVu = 165 ib 5 GeV.Q A common approach is to use the 
minimization conditions and RG equations to solve for the inputs, given a phenomenologically 
allowed set of weak-scale outputs. In the case of a high messenger scale, however, no easily 
invertible solution for the RG equations is available. Instead, we simply choose to tackle the 
problem numerically. After running down the RG equations and minimizing the Higgs potential 
once, we iterate this procedure, each time adjusting the value of the parameter B to fix the overall 
scale of the VEVs and masses, while simultaneously changing the dimensionless couplings to 
correctly reproduce the top quark mass. This iteration process, in fact, converges fairly quickly. 

Using the procedure above, we perform a numerical scan of a large portion of the parameter 
space. We study the low-energy particle spectrum for various messenger scales A, numbers of 
messengers n, and values of the couplings A and k. It is interesting to note that it is very easy 
to generate non-zero VEVs for H^, H^, and A^, even when is a small positive number. This 
is because the terms \XvdVu — kx^\'^ and AxXvdVuX, when XvdVu and k are of the same sign, both 
"push" the VEV of the real component of the singlet away from the origin. Unfortunately, 
we find that, for any choice of values of the input parameters, there are always particles with 
unacceptably small masses. To illustrate the situation, we present in Table |I| our numerical 
results for several representative points in different "corners" of the parameter space. The first 
two points represent the typical situation for the case of the low-energy GMSB, the next two 
are representative of the case of the high-energy GMSB, and the last one explores the extreme 
case of A = 10^^ GeV. Points 1 and 3 have relatively large values of k, while points 2 and 4 
have /c -C 1. Notice that in the table we did not consider a similar limiting case for A. This is 
not a coincidence. It turns out that, for A ^ 0.2, the dominant term in the potential is Vn, and 
tan/3 is forced to values very close to one. In this case, in order to correctly reproduce the top 
quark mass, one is forced to choose hf at the weak scale such that ht hits the Landau pole below 
the GUT scale. We have chosen to list only the cases where the couplings in the superpotential 
remain perturbative up to the GUT scale. We make, however, no such assumption in our analysis 
in the next subsection. 



"Notice that this number is not equal to the top quark pole mass, the experimentally measured quantity, 



because of QCD corrections. The relationship between the two is given, at 1-loop, by nipoie = m{l 
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Table 1: The numerically determined NMSSM parameters for five sample points in the parameter 
space. Here m/j. and refer to the eigenvalues of the scalar and pseudoscalar Higgs mass 
matrices respectively, and mg denotes the mass of the right-handed selectron. The values of A, 
k, and hf are given at the weak scale. All the other quantities have been defined earlier in the 
text. 





Input Parameters 


point 


A (GeV) 


A 


k 


B (TeV) 


n 


ht 


1 


5 X 10^ 


0.25 


0.1 


6.4 


1 


1.12 


2 


10^ 


0.28 


3 X 10-^ 


3.6 


3 


1.08 


3 


1012 


0.32 


0.3 


0.99 


10 


1.07 


4 


10^2 


0.25 


3 X 10-4 


6.0 


1 


1.11 


5 


lO^'' 


0.28 


0.3 


6.9 


1 


1.07 





Soft SUSY-breaking Parameters at the Weak Scale 


point 


(GeV^) 


(GeV^) 


(GeV^) 


A, (GeV) 


Ak (GeV) 


1 


-2.4 X 10^ 


5.3 X 10^ 


4.6 


-1.5 


-4.0 X 10-=^ 


2 


-2.8 X 10^ 


5.7 X 10^ 


6.8 


-2.6 


-6.2 X 10-=^ 


3 


-3.1 X 10^ 


4.8 X 10^ 


29 


-11.4 


-0.15 


4 


-2.5 X 10^ 


6.8 X 10^ 


12 


-8.0 


-0.11 


5 


-2.9 X 10^ 


1.0 X 10^ 


-8.1 


-9.4 


-6.0 X 10-3 






Field VEVs 


Particle Masses 


point 


tan/3 


X (GeV) 


Ms (GeV) 


me (GeV) 


rrih, (GeV) 


rriA, (GeV) 


1 


1.59 


-3.7 


61 


32 


85, 39, 35 


51, 1.8 


2 


1.84 


-3.7 


103 


35 


87, 48, 38 


48, 0.2 


3 


1.97 


-40 


94 


36 


87, 53, 28 


76, 25 


4 


1.63 


-14 


57 


34 


85, 43, 37 


44, 0.5 


5 


1.88 


-49 


66 


40 


88, 50, 27 


71, 24 
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It can be seen that, in all the cases presented in Table |l|, there are particles with unacceptably 
small masses. The result for the low-energy GMSB is not new and has been known for several 
years 0. On the other hand, the situation with a high messenger scale had not been quanti- 
tatively studied in the literature to date. One expected feature that we indeed see in points 3 
and 5 is the increase of the pseudoscalar Higgs boson mass with A. This happens because the 
magnitude of Ax, generated by running the RG equations, increases with the messenger scale, 
and it is Ax that breaks the U{1) symmetry of the potential, as discussed before. Another result 
that could have been anticipated is the smallness of the mass of the light pseudoscalar Higgs 
when k <^ 1 (points 2 and 4). This is due to the Peccei-Quinn symmetry, which is restored 
in this limit. What is surprising is that raising the messenger scale by 10 orders of magnitude 
does not bring any other significant changes to the particle spectrum. The masses of the gluino, 
right-handed selectron, and scalar Higgs boson still remain small. 



4.3 Analytical Considerations 

In this subsection we present a rather simple semi-analytical argument which explains why there 
can be no phenomenologically acceptable solution to the NMSSM with the GMSB. We show 
that if one assumes that such a solution exists, one arrives at a contradiction. We also explain 
some of the features of the numerical solutions presented in the previous subsection. 

Suppose that for some point in the parameter space an acceptable solution exists. The 
problem to be addressed is the smallness of the selectron, gluino, and Higgs masses. We choose 
to base our analysis on the right-handed selectron mass constraint. The magnitude of rrig is 
directly proportional to the size of the 5-parameter. In our numerical procedure the value of B 
is chosen in such a way that v = 174 GeV. A typical value of B obtained in this way yields a very 
small selectron mass (mg ~ 35 GeV), gluino mass (M3 ^ 100 GeV), and soft SUSY-breaking 
masses for the Higgs bosons {fn'jj^ ~ —3000 GeV^, m|^^ ~ 500 GeV^). 

It is, therefore, obvious that the only chance of obtaining an acceptable value for the selectron 
mass is to raise B, by a factor of three or more, and try to arrange the other parameters in 
such a way that \Jv'^ + remains 174 GeV. Since B feeds into all soft SUSY-breaking masses, 
their absolute values will also increase. For example, imposing mg > 80 GeV forces fn'jj^ < 
— (215 GeV)^ for a messenger scale of 10^^ GeV. For different messenger scales the bound becomes 
even more stringent, as shown in Fig @(a). 

To determine the consequences of raising the soft SUSY-breaking masses, we analyze the 
Higgs potential (Eq. ([4.5| )). The extremization conditions at tree level are 

kx^)Xvu + 2X^x\d + 2ml^Vd - 2AxXvuX + ^ ^ 2vd{vl - vl), (4.8) 

0'2 I 2 

kx^)\vd + 2AW„ + 2m]j Vu - 2Ax\vdX + ^—-^2vu{vl - vl), (4.9) 
kx^){-2kx) + 2xX'^{vl + vl) + 2m%x - 2AxXvdVu - 2kAkX^ . (4.10) 

The first two equations (Eqs. (|4.8|) and ([4.9|) ) closely resemble the corresponding ones in the 
MSSM case. In fact, the only difference in the NMSSM is the presence of the first term on the 
right-hand side of Eq. ( [4. 81) and Eq. (|]9|). This term originates from = \XvdVu — fcx^p and 

is, therefore, absent in the MSSM. Dividing Eq. ( ^^ ) by Eq. ( [4.9| ) by Vd, and subtracting the 
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Figure 2: (a) Lower bounds on \mffj'^^'^ and \mjjj^^'^ as a function of the messenger scale A from 

the selectron mass constraint mg > 80 GeV. Here n = 1, ht = 1.07, k = 0.3 and A = 0.29 at 
the weak scale. These bounds do not change for different values of k or A. The other plots show 
typical values of (b) Ax, (c) Ak, and (d) m^, for the same choice of parameters that yielded (a). 
The values of these parameters do not change significantly for different values of k or A. 
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two expressions, we can cancel out this term. As a result, we obtain 



2^2 



Ml 



m 



H, 



+ m 



m 



Ml 



cos 2/5 



+ 



m 



tan^ (3-1 



(4.11) 



Note that this equation is identical to Eq. ( p.5|) with /i = Ax. To obtain the NMSSM analog 
of Eq. (|3.6|) we divide Eq. (|4.^ ) by Vd, Eq. ( [4.9|) by Vu and add the two. Solving for A\\x 
{A\\x = m\) we find: 



AxXx = [m 



, ,9 9, sin 2/5 
+ m^^+2AV)^-^ 



+ A(At;dt;„ - kx"^) 



(4.12) 



Eqs. ( ^.61 ) and ( ^.12| ) differ only by the contribution from 

Eq. ( [4.11|) states that the value of the effective //-parameter generated in this model is subject 

M|/2, which, if one imposes m\j < —(212 GeV)^, 



-m 



to a rather stringent bound: A x > 
translates into Ax > 200 GeV. Notice that the origin of this bound is the same as of the bound 
on the size of the //-parameter derived in Section ^, since the condition given by Eq. dill 



IS 



the same in both cases. In the present case, however, the bound is stronger because tan [3 is no 
longer a free parameter but is determined by minimizing the Higgs potential. 

So far we have only looked at the first two extremization conditions. We now turn our 
attention to Eq. ( [4.10| ). Solving for x^ in Eq. (|4.11| ), one can rewrite Eq. (|4.10| ) as 




Hu 



ia.v? (3 



Mr 



tair (3 



\v'^{ksin2(3 -\) 



m 



2 , . X 2sin2/? 
^ + AxXv 



2x 



+ kAkX. (4.13) 



While we have shown that phenomenology requires the expression in parenthesis on the left-hand 
side to be larger than (200 GeV)^, the terms on the right-hand side are all much smaller, because 
m%, Ax, and Ak are zero at the messenger scale and the effects of the RG running are relatively 
small (see Fig. ^). This means that the above equation can never be satisfied unless k <^ X. 

An immediate consequence of the k ^ limit is that the mass of the lightest pseudoscalar 
Higgs goes to zero, as it becomes a Nambu-Goldstone boson. (It is for this reason that k was 
introduced in the first place.) Furthermore in the limit of large /i and small k the determinant of 
the scalar Higgs mass-squared matrix becomes negative, which means that the extremum point 
given by Eqs. ( ^4.8| - [4.10|) ceases to be a minimum. To show this we first derive a relationship 
between A; and sin 2/?. That relationship can be derived from Eq. ([4.11|) and Eq. ([4.12|) . Neglecting 
M|, AxXx, and A^f^ in comparison to fn'jj^ and fn'jj^, we find that 



sm2/3=y ^ 

^l + r^d) 



1 + r 



V 




(4.14) 



Here r = —{rri^j^ + fnjjj / {m'jj^ — mjjj^ For k <^ X, Eq. {\i.l4\) reduces to 



sin 2/? ^ 2- 
A 



-m 



Hu 



^Ha-^Hu. 



(4.15) 



ttin deriving Eq. ( 4.14 ) it was necessary to assume that sin 2/3 > k/X. This translates into two requirements: 
r > and k/\ < 1. We conclude that for large soft SUSY-breaking Higgs masses-squared it is necessary to have 
\>k. 
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Equipped with the last result, we consider the determinant of the scalar Higgs mass-squared 



matrix. The full expression for it is given in [Appendix B| ; here we only need to identify the 



leading terms. We are interested in the case n > Xv, and, as we have argued before, the soft 
trilinear couplings A\ and Ak can be neglected. For this reason, the dominant terms will be the 
ones containing the highest power of fi: 

det Ml^i,, ^ Ai^§^ (-4^^' + ^'9' + cosmf + Sk'X' sin(2/5)) . (4.16) 

That these are indeed the largest terms was checked numerically. 

Taking into account the fact that k and sin 2/3 are proportional to each other for small k, 
one can easily see that, in the limit A; ^ 0, the first term dominates and the determinant is 
neg ative.0 

This completes our argument, and we are now able to state that there can be no phenomeno- 
logically viable solution in the context of the NMSSM. We could have also based our argument 
on the gluino mass bound. The experimental constraint M3 > 190 GeV translates into the 
requirement m|^^ < —(212 GeV)^ (assuming n = 1), and the rest of the argument follows un- 
changed. Notice, however, that the bound on m'jj weakens if the number of messenger fields is 
taken to be very large. 

We now turn to the issue of interpreting the numerical results of the previous subsection. 
We would like to understand, for instance, why the values of the singlet VEV x in Table are 
always smaller than the VEVs of the Higgs doublets and, furthermore, why x is only several 
GeV for a low messenger scale. 

The answer comes from considering the extremization condition for x: 

2Px3 + A (A - A; sin(2/5)) v^x - ^^^(^Z^)^^ ^ q , (4.17) 

where we omitted the terms and kA^^x (|m^| <^ A^f^ for all the points in the table). For 
most of the parameter space the cubic term in x can also be neglected, giving 

Ax sin(2/j) 



2(l-fsin(2/5) 

Thus the smallness of x is related to the fact that A\ is small. The above approximation holds 
as long as 

A2 2 fl - ^ sin 2/5)^ 
Al<\'v% ^ / , (4.19) 

k^ sm Ip 

which is not satisfied only for point 5 in Table ^ For point 5 the value of x can be approximated 
by 



2 A sin 2/5 



1/3 



Again x < v and therefore Ax ^ 175 GeV. 



^Because sin 2/3 (xk/\ there is no ambiguity with sign redefinitions of A or fc in Eq. (4.16) 



20 



Knowing that x is small in this model we can derive another interesting relation. Neglecting 
all the terms containing x in Eq. (|4.12|) , we obtain: 



AV^-(m^^+m^J. (4.21) 

This explains why the values of the soft SUSY-breaking masses for the Higgs bosons are so 
similar for very different values of the messenger scale. 

Finally, we can say a few words about the scalar Higgs boson masses. In the limit of small 
X (and hence small /i), the dominant term in the determinant of the scalar Higgs mass-squared 
matrix (see [Appendix is 



(4.22) 



Taking into account the fact that, for small x, /x = Ax ~ (see Eq. ([4.18|) ), the equation above 
gives: 



32sin(2/3) 



,2 _2 _2 



(4.23) 



This explains why changes in the messenger scale have almost no effect on the product of the 
scalar Higgs boson masses (see Table |I]), as long as A is unchanged. For sin(2/5) ~ 0.8 — 0.9, 
which is what we typically find in this case, Eq. ( [4.23| ) gives a "geometrical average" value of 



the scalar Higgs boson mass of only about 50 GeV. This means that, as long as x is small, the 
model necessarily yields phenomenologically unacceptable Higgs boson masses. 



5 Possible Modifications to the NMSSM 

In this section we reexamine the expressions derived in Section ^ and attempt to modify the 
NMSSM to make it phenomenologically viable. We study several possibilities and comment on 
the problems that arise. Overall, we find none of these possibilities entirely satisfactory. 



5.1 Extra Vector-like Quarks 

We want to modify the NMSSM in a way that allows one to avoid the conclusions of Section ^. 
Recall that the crucial step in our analysis there was the observation that Eq. ( [4.13| ) could 
not be satisfied: the left-hand side was always greater than the right-hand side. To obtain a 
consistent solution one has to somehow make both sides equal. One possibility is to make 
of the same order of magnitude (and sign) as . That could be accomplished by coupling the 
singlet to some new fields and arranging the parameters in such a way that the SUSY-breaking 
mass-squared of the singlet is driven sufficiently negative. This idea was first proposed by Dine 
and Nelson in Ref. 0, who introduced new color-triplet fields g' and q' and coupled them to A^. 
The corresponding superpotential is 

W = KQHuu' + hdQHd" + KLHae'' + XNHdHu -hn'' + X.Nq'q' . (5.1) 

According to Eq. ( |2.6| ), the scalar components of q' and q' acquire large SUSY-breaking masses, 
which can drive mj^ sufficiently negative. 
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Agashe and Graesser in Ref. did a quantitative study of this scenario for the case of 
the low-energy GMSB. They showed that it is indeed possible to generate a large negative 
m%, in the range —(150 GeV)^ to —(200 GeV)^, and further demonstrated that, with of 
this magnitude, one can choose the input parameters in such a way that v = 174 GeV and 
all experimental constraints are evaded. They also pointed out that in this scenario the input 
parameters need to be fine-tuned in order to reproduce the above value of f . In what follows we 
give a set of input parameters that yields an acceptable particle spectrum, and then proceed to 
analyze the sensitivity of the Higgs boson VEVs to the NMSSM coupling constants. We clarify 
the origin of this sensitivity and also extend the analysis to the case of the high-energy GMSB. 

As an example of an allowed solution, we consider the case of the low-energy GMSB with 
B = 50 TeV, n = 1, and A = 100 TeV. For m% = -(190 GeV) 2, to correctly reproduce 
Mz and rrit we take ht=0.99, k = —0.045 and A=0.11 at the weak scale. We find that tan/3 
equals —2.9 for this point. Because the magnitude of the product Bn is now quite large, the 
masses of the gluino and right-handed selectron are safe: M3 = 477 GeV, mg =93 GeV. The 
vacuum expectation value of the singlet is also large, a:=2.97 TeV, which, as was argued earher, 
is required by Eq. ( [4.111 ). The eigenvalues of the scalar Higgs mass matrix are 404, 270, and 90 



GeV, and those of the pseudoscalar Higgs mass matrix are 400 and 6.7 GeV. The last number 
appears alarmingly small at first sight but, as shown in Ref. P2|, has not been excluded. The 



reason is that the corresponding eigenstate a is almost a pure singlet: 

\a) = Omi\Hd) - OMl\Hu) - 0.999|A^) (5.2) 
The quantitative criterion given in Ref. |^ , based on the constraint from the T —>■ a'-f decay, is 

sin 2/5 tan/? 



(250 Gev) 



^ + sin^ 2(3 



< 0.43, (5.3) 



and for the parameter set above the left-hand side equals 0.15. 

In this scheme it is, therefore, possible to find a point in the parameter space which leads to 
a phenomenologically viable solution. Unfortunately, as we already mentioned, this solution is 
very sensitive to the choice of the superpotential coupling constants A and k. In the remainder 
of this subsection we discuss this issue in detail. 

The values of the parameters for the set that we have just described had to be chosen in such 
a way that the top quark and Z-boson masses were fixed at their known experimental values. 
It is interesting to investigate what values of Mz would be predicted for a generic choice of the 



parameters. In Fig. ^ we plot the magnitude of the quantity v = yv'^ + as a. function of A 
and k. The figure shows that small changes in both A and k lead to large changes in v. This is 
very similar to the situation in the MSSM which was considered in Section There we showed 
that the value of the /i-parameter had to be chosen very carefully in order to yield the correct 
value of V. In the present case, the points in the parameter space that correspond to values of v 
around 174 GeV lie in a very thin band on the X — k plane. Also notice that, for this range of A 
and k, the slope is the steepest. (See [Appendix C| for comments on this point.) 



It is possible to perform the same type of analysis for a higher messenger scale. The same 
problem is found in that case as well. In Fig. ^ we plot the dependence of v on A for fixed values 
of k. For comparison, the curve for A = 10^^ GeV is plotted next to the curve for A = 100 TeV. 
From the slopes of these curves one can determine the degree of sensitivity with respect to A, 
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Figure 3: The value of f = yf J + 1;^ as a function of A and k. The inputs are n = 1, = 
-(190 GeVf, B = 50 TeV, A = 100 TeV, ht = 0.99. 



' — 1200 p — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r 




Figure 4: The dependence of v on the value of A for the high- and low-energy GMSB. The other 
input parameters are the same as in Fig. ^ 
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using the definition in Section ^ The degree of sensitivity, given by d(logf )/d(log A), is 2% for 
the low-energy curve and 1% for the high-energy curve. Our numerical results agree with those 



in Ref. for the low-energy GMSB if the same inputs parameters are used. 

In order to understand this behavior, we once again turn to the extremization conditions 
Eqs. ( |4.8[ - [4.10| ). First, we present some qualitative observations. Recall that phenomenology 



requires \x\ to be rather large (of the order ^J\mjjJ/\ ^ 1 TeV), while v has to remain "small" 
{v = 174 GeV) to correctly reproduce Mz- As a result, the terms containing high powers of x 
and the terms containing m? [i = Hd, Hu, N) dominate, while the terms with Vu and Vd are not 
fixed, and have to absorb the residual difference between the dominant terms. Therefore, small 
percentile changes in the dominant terms can result in large percentile changes in the Higgs 
boson VEVs. This is to be contrasted with the situation in the previous section, where A^f ^ was 
tied to the value of the sum mj^^ + m'jj^ (see Eq. ( [4.21|) ). 



Next, we try to identify the main source of this sensitivity. We first consider the dependence 
of w on A for fixed B, k, and hf. One can use Eq. ( [4.11| ) to solve for v"^ and then isolate the 
largest contribution to dv/dX. 



dv^ 4 



dX 



2 o\2 ^ I mjj^- mjj^ tan^ p\ d tan (3 



-2Xx' - 2X^x— + 



dX dtanP \ tan2/5-l / ^A 
1 dm\j^ tan^/5 <9m^^ 
tan2/5-l dX " tan2/3-l OX 



(5.4) 



Using the data that led to Fig. |^, we numerically evaluate the derivative around the point 
A = 0.11, k = —0.045. The following are the results of evaluating each of the terms on the 
right-hand side, respectively: -1.4 x 10^ -2.4 x 10^ -3.9 x 10*^, -1.2 x 10^ 1.9 x 10^ (GeV^). 
The largest term is the first one, the next two terms combined provide a 45% correction, and 
the derivatives of the soft SUSY-breaking masses can be completely neglected. In [Appendix D| 



we show how these numbers can be understood by studying the minimization conditions. 

The fact that the dominant contribution to dv/dX comes from the first term in Eq. 
has a very important implication. It means that the problems of cancellation in the NMSSM 
and the MSSM are not merely similar, but have exactly the same origin. Indeed, Eq. ( [4. Ill) 



is the same as Eq. ( p.5|) , and, because in the NMSSM v depends on A mainly through the 
combination Ax, which plays the role of the /x-term, the two models require roughly the same 
degree of cancellation. The degree of cancellation quoted in Section ^ for the MSSM is most 
conservatively 16%, but this is so because one can choose tan /3 freely in the MSSM. On the other 
hand, tan (3 is determined by minimizing the potential for the NMSSM and cannot be chosen 
arbitrarily to ease the cancellation. For the value of tan/3 which we obtained in the NMSSM, 
the degree of cancellation is actually comparable (order a few percent) in the MSSM. The small 
difference between the two models is due to the dependence of x and tan P on A. 

We have discussed the A dependence of the Higgs boson VEV, and now turn to the k de- 
pendence. Fig. |] shows that the points that yield v = 174 GeV form an almost straight line on 
the X — k plane. It can be shown (see [Appendix L)| ) that in order to keep v constant one has to 
change k and A according to Ak/k = A A/ A. The sensitivity of f to is, thus, related to the 
sensitivity of v to A, which, in turn, originates from the need to carefully choose the /i-parameter 
in the MSSM as discussed in Section 

To summarize, we have shown that this model requires a very particular choice of parameters 
to yield the correct Z-boson mass. Furthermore, we explained that the sensitivity of the Z-boson 
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mass to the NMSSM couplings has the same origin as the sensitivity of the Z-boson mass to the 
value of the /x-parameter in the MSSM. We emphasize that the problem is present for both high 
and low messenger scales, simply because the bound on the /^-parameter does not weaken as one 
raises the messenger scale. 



5.2 Hypercharge D-term 

Next, we investigate what happens if the D-term contributions described at the end of Section]^ 
are included. First, we consider the case of the NMSSM with no extra particles added. We try 
to determine if, by introducing the D-terms, it is possible to make v smaller. If that happened, 
V could be rescaled back by increasing B, and that would raise all masses in the model, as 
desired. We find that this is not the case. Upon adding the D-terms both tan j3 and (x) change, 
but f J + vl, curiously enough, remains virtually constant. This happens because, in the limit 



X 



^ f ^, f ^ is constrained by Eq. ( |4.21| ), and the change mjj^ ''^Hd ~ ''^Hu ~^ "^Hi, + 
preserves the quantity m'jj + . 

The next question to ask is whether the D-terms can decrease the degree of cancellation for 
the case with q' and q' added. The answer is again negative and the reason can be seen from 
Eq. ( |4.12| ). Recall that the degree of cancellation is controlled by the magnitude oix^. As long 



as AxXx and X^VdVu can be neglected compared to rn?^^ + rn?Hu^ ( |4-12|) yields 



x'^-'-^ (5.5) 



and the relevant quantity is again rn^j + m?^ 



5.3 Large Trilinear Couplings 



They consider a 
A <^ 1 and the value of Ax is large at the 



At last, we consider the scenario proposed by Ciafaloni and Pomarol 
modified version of the NMSSM, where = 0, 
messenger scale. Their model also contains, in the potential at the weak-scale, a linear term 
in N which is generated by tadpole diagrams and solves the problem of a light pseudoscalar. 
They find that the requirement of the positivity of the determinant of the scalar Higgs boson 
mass-squared matrix is very restrictive. We repeat part of their analysis to determine if their 
choice of parameters could indeed lead to a phenomenologically viable electroweak symmetry 
breaking spectrum. Note that, as far as the following is concerned, their model is identical to 
the NMSSM. 

The full expression for the determinant can be found in [Appendix B| . In the limit of ^ 
and A — > 



scalar 



A|t;^M|A^sin'(4/j) 



l + y- 



cos 



'3. 

Ml 



+ 1 



y 



cos2(2/3) 



(5.6) 



where we introduced a variable y = g'^mj^ / {2M^ 
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^/J to conform to the notation used in Ref. 
From the extremization conditions for the potential, Eqs. (4.S- 4.ld| ), one can show that 
sin(2/3)/(2(l -|- y)). There are two intervals of y over which the determinant is greater 



than zero. One interval is where both the expression in the brackets and the denominator are 
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positive. It is given approximately by the following bound on \y\: 

'1/2 



\y\ < 



( 

cos 2/5 1 + ^ 



Ml 



■^1 



The other interval, not mentioned in |31], is approximately 

denominator and the bracketed expression are negative. 

The first interval, for Ax > Mz, corresponds to rather small values of and 

/i ~ ^Ax sin(2/5) . 



(5.7) 



where both the 



Using this equation together with Eq. ( |4.12| ), one can derive the following result: 



(5.9) 



The above equation is impossible to satisfy in models with the GMSB, because the combination 
(m^^ + trijj^) is always negative at the weak scale for the messenger-scale boundary conditions 
given by Eqs. ( |2.5| ) and ( |2.6| ). To satisfy Eq. ( p. 91 ), a drastic modification of the boundary 
conditions would be required. 

We now turn our attention to the second possibility. It requires a relatively large negative 
value of the singlet soft SUSY-breaking mass-squared: < —2/{g'^) x M^ = —(132 GeV)^. 
This value is impossible to generate unless, as before, one introduces fields q' and q' and couples 
them to A^. Even with the introduction of these fields, if = 0, A -C 1, the extremization 
conditions cannot be simultaneously satisfied. This can be seen in the following way. For k = 
Eq. ( [4.13 ) takes on the form 

AxXv^ sin(2/?) 



X 



m 



(5.10) 



which implies x 
as A — > 0. 



as A ^ 0. This is incompatible with Eq. ( [4.11[ ), which requires that x oo 



6 Conclusion 

We studied the issue of electroweak symmetry breaking in models with the gauge mediation 
of supersymmetry breaking (GMSB). We first reviewed various proposals in the literature to 
generate the /x-parameter of the MSSM with the same order of magnitude as the soft SUSY- 
breaking parameters such as squark, slepton, and gaugino masses. We find that most of them 
require small parameters which are accidentally of the same magnitude as the loop factors, 
cancellation of the kinetic mixing terms at the level of 10~^, omission of interactions allowed by 
symmetries, or many new degrees of freedom not motivated otherwise. 

Even if one could generate the /x-parameter with the same order of magnitude as the soft 
SUSY-breaking parameters, it has to have particular values to reproduce Mz = 91 GeV. We 
studied this question numerically and found the following. The current experimental lower 
bounds on superparticle masses limit the overall scale of SUSY breaking from below, which in 
turn limits m|^^ < from above [i.e., {rnj^J from below). To reproduce Mz, fi"^ needs to cancel 



26 



(too-negative) rn'jj^ and is hence bounded from below. Therefore, there is some cancellation 
required between /i^ and m'jj^. Even with the most conservative set of parameters, we found 
that a cancellation of 16% is necessary. The situation is worse for most of the parameter space. 
This situation was contrasted to the supergravity scenario where the current experimental lower 
bounds on superparticle masses do not require a significant cancellation among parameters. 

The simplest mechanism to generate the //-parameter would be the NMSSM, the minimal 
extension of the MSSM without dimensionful parameters in the superpotential. The NMSSM 
is known not to work with the low-energy GMSB, but there was hope that it may work with 
higher messenger scales. We have shown that this is unfortunately not the case. The current 
bounds on the superparticles masses are already strong enough to exclude the model completely. 
We presented a semi-analytic discussion to clarify why the NMSSM fails. 

We also discussed various possible modifications to the NMSSM and whether they could lead 
to a viable electroweak symmetry breaking. The introduction of extra vector-like quarks coupled 
to the NMSSM singlet produces a large negative mass-squared for the singlet, and leads to a 
viable electroweak symmetry breaking. One needs to adjust the parameters to a few percent, 
which is comparable to the MSSM case for the same tan/? range. A Fayet-Illiopoulos D-teim 
for U{1)y does not improve the situation. 

The overall prospect of electroweak symmetry breaking with the GMSB remains unclear. We 
hope our detailed investigation prompts further studies on this issue. 
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Appendix A The Renormalization Group Equations of 

the NMSSM 

In this appendix we list all of the RG equations for the NMSSM, at 1-loop [^. These are the 



equations used, in section 4, to determine the coupling constants and SUSY-breaking parameters 
of the NMSSM at the weak scale, given their values at the messenger scale. 

IGTT'j^g' = llg'^ (A.ll) 

167r^^^2 = gl (A.12) 

167r2^^73 = (-3)^/1, (A.13) 

Wn'j^h, = {Qh',+hl + X'-fg"-3gl-fgl)h„ {AM) 
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rl+ 

ox 


= {6hl + hi + hl + X'- \f - -igl - 

y 


O 


(A.15) 


at 






(A.16) 


dt 


= (4A' + 2A;' + 3/i^ + 3/i^ + /i^-^''- 




(A.17) 


dt 






(A.18) 



In the above equations g' is the U{1)y gauge couphng; exphcitly g' — e/ cos 9w- g2 and gs 
are, respectively, the weak and strong couphng constants. One defines gi to be the hypercharge 
couphng constant in the GUT normahzation, i.e. gi = ^f^g' and ai = |q;'. Gauge couphngs at 

the messenger scale are defined in such a way that they match their experimental values at the 
Z-mass. We only consider the effect of third generation Yukawa couplings, namely, ht, hb and 



dt 



4{^g''M, + ^-glM, + ^-glM,), (A.19) 



levr^^Ad, = &hl{l + 5aMh + '^hl5a3At + 2hl5azAr + 2\^Ax 



dt 



i{^g"M, + ^glM2 + ^^g'^s), (A.20) 



IQn^^Ae^ = 2hlil + 35a3)Ar + 6hlA, + 2X^Ax 



dt 



- Q{g'^Mr+glM2), (A.21) 
IGtt^-^^a = S\^Ax-Ak'^Ak + QhlAt + QhlAj, + 2hlAr 

- 2{g'^Mi + 3glM2), (A.22) 
IGTr^-^Afe = 12{k^Ak-X^Ax). (A.23) 

Ai are the soft SUSY-breaking trilinear couplings, given in Sections 2 and 4. Note that we only 
consider third generation trilinear couplings, namely Afht = A"^, Abhi, = A^^, Arhr = Af^. Mi 
(i=l,2,3) are the soft SUSY-breaking gaugino masses and they evolve, at one loop, identically 
to ai. Explicitly 

Ml gi • ^ ' 

where gx is the value of all gi at the GUT scale, while Mi is the common gaugino mass at the 
GUT scale. 

167r';^m|^ = '^^ash^iml^ + m^^ + ml + A^) + 2Sa3hl{ml^ + mj,^ + ml + Al) 



- Si^g'^M! + \glMi + \glMl) + \g'% (A.25) 
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8{^g'V, + ^gM) - ^g'X (A.26) 



dt 

- 8{^g'V,+^gM) + yX (A.27) 



- 8{\g'V,+^glM',)-g% (A.28) 

levr'^m^^ = "^Sashliml^ + ml^+ml + Al) 

- 8(?'2M2 + 2^7'^e, (A.29) 



+ 2A2(m^^ + ml^ + m% + Al) - 8i}-g'''Ml + \glMl) - ^7'^^, (A.30) 

167r^^m|„ = 6/i?(m|^ + m^^+m? + A2) + 2A2(m^^+m^^+m^ + A2) 

- 8i}-^g'^Ml^\glMl)^g'% (A.31) 

levr^^m^ = 4A2(m^^+m^^+m^ + A2)+4A;2(3m^ + A2). (A.32) 

^ is the hypercharge-weighted sum of all soft SUSY-breaking masses- squared 

i = Y^Y,m\, (A.33) 

i 

where i runs over all scalar particles. With the boundary conditions in Eqs. ( p.5| , |2.6|) , ^ = and 
remains zero throughout the RG evolution. All soft SUSY-breaking mass-squared terms were 
taken to be diagonal. Again, we only consider the running of third generation soft SUSY-breaking 
masses- squared, is defined in Section ^. 



Appendix B Scalar Higgs Mass-Squared Matrix 

In this appendix we explicitly show the 3x3 scalar Higgs mass-squared matrix of the NMSSM. 

^scalar 2 OViOVj 2 

( fvl + (Aa + ^) ^ (4A2 - f)viV2 - 2/i (A^ + ^) 4A/ii;i - 2A^\v2 - AkfiV2 \ 
(4A2 - g^)v,V2 - 2fi + f) fvl + (a, + ^) ^ -2A^\v^ - Akfiv, + AXfxv^ 



(B.l 
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where Vi for i = 1,2,3 corresponds, respectively, to v^, Vu and x. All other parameters were 
defined in previous sections. 

The determinant of the matrix above can be evaluated explicitly, and its full expression is 
given bellow. Various limits of this determinant are considered in the body of the paper. 

det A<L/ar = 

SzA'^/i sm(2/y) ^ 

128kv^X^fi^ - 32Akk^X^fi^ - 128AaA;=^AV - 256AxX^fi^ - 256A;AV^ + SAxv^X^ cos(4/5) + 
32A^AV^ cos(4/5) + GAAxv^X'^fi^ cos(4/5) + 32AkAxkX^fx^ cos(4/5) + IGOA^fcAV^ cos(4/3) + 
128fcv2A V cos(4/5) + 32AfcA;2AVcos(4/5) + 128Ax A;^A V cos(4/3) - 2AAt;^A^ cos(8/?) + 
SAxv^X'^f + 32Axv^X^fi^f - IGAkAxkX^fi^f + QAkv^X^i?f - IGAkk^XfiY + 
GAAxk'^Xjj.Y + Q^k^lJ^^f - 4Axv^X'^cos{A(3)f - 32AAf^AV cos{4(3)f - 
16AkAxkX'^fx^cos{AP)f - GAkv^X^fx^ cos{AP)f - IGA^k'^Xfj,'^ cos{AP)f + 
64AAfc^A/ cos{4:p)f + Uk^ii^ cos{A(3)f + Axv^X'^ cos{8(3)f + ASAlv^X'^fx sin(2/5) + 
24AfcA;i;2AV^sin(2/3) + UOAxkv^X'^fx^ sm{2p) + 256^^ A V sm(2/5) + 96t;2A V sm(2/3) + 
768AaA;A V sin(2/5) + 512PaV sin(2/3) - 12A2t;2A5^^2 ^^^^2/?) - 12Afcfct;2AV^^^ sin(2/3) - 
60Axkv^X'^fj.^fsm{2P) - 16k\^X^fi^f sm{2P) - 48t;2A sin(2/5) - WAlv^X'^ fxsm{6p) - 
SAkkv'^X^l? sin(6/3) - AQAxkv'^X^ ^i^ sin(6/5) - 32v'^X^ii^ sin(6/3) + 4A2t;2A5^^2 gin(^Q^) + 

4AfcA;t;2AV^^sin(6^) +4AAA;t;2AV^^sin(6/?) - 16A;2t;2A sin(6/5) + IQv^X^ ii'^f . 

(B.2) 

All parameters were defined previously. Recall that /i = Xx. 

Appendix C Comments on Naturalness 



We have studied the NMSSM with extra vector-like quarks in Section [ST^ and discussed that 
the model requires a delicate cancellation among independent parameters. In this appendix, we 
make further comments on the naturalness of this model. 

From Fig. ^, one can easily note that not only does the experimentally allowed value of v lie 
on a steep region of the parameter space, which requires a degree of cancellation of order 1%, 
but it lies on the steepest region of the parameter space. 

One may, therefore, try to address the following question: if all parameters are kept fixed (and 
this choice of parameters yields an experimentally allowed spectrum) except one {e.g. A), what 
is the likelihood of obtaining a certain value of v upon a random choice of the free parameter? 



In other words, what is the probability P{v) dv of finding the value of yf^ + v1 between v and 
v + dv given a random choice of A? This line of reasoning is related to the definition of fine-tuning 



introduced by Anderson and Castano It is easy to note that 



pwoc(g\ (C.I) 

This "probability density" is plotted in Fig. ^. Note that we restrict A to lie on a range where the 
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1000 



Figure 5: The probability densities of finding specific values of v in the NMSSM with extra 
vector-like quarks upon random choices of A. All other parameters are the same as in Fig. ^. 
The probability densities are normalized so that P{v = 174 GeV) = 1. 



same "qualitative" physics is obtained, that is, electroweak symmetry is broken and tan(/3) > 1. 
The plot has been normalized in such a way that P{v = 174 GeV) = 1. 

We note that, in some sense, the probability of living in our universe is smaller, if this model 
is to be taken seriously, than the probability of living in a universe where v ~ 600 GeV by a 
factor of three. One can turn this picture around and say that the NMSSM, with the above 
choice of parameters, "prefers" (or predicts) v ~ 600 GeV. 

This does not happen in the MSSM. The analog of Fig. |^ would be Eq. ( ^.51 ), which is 
a straight line (M| = M|(/i^)) if all parameters except /x^ are kept fixed. In the language 
introduced above, the MSSM does not "prefer" (or predict) any particular value of M|, that is, 
the "probability density" of M| upon random choices of /x^ is fiat. 



Appendix D The Dependence of the Higgs VEVs on the 

couplings of the modified NMSSM 



In Subsection |5.1| we showed that the values of the Higgs boson VEVs were extremely sensitive 
to small variations of the superpotential couplings A and k. These variations were evaluated 
numerically after Eq. (|5.4|) for one particular set of A and k. In this appendix we study this 



issue analytically and show how one can estimate the effects of small variations AA and Ak on 
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V, tan (3, and x. 

We will use the following three equations, derived in Section 



4 



-AV + /, (D.l) 



2A;V ~ -mj^ + X{ksm2l3- X)v\ (D.2) 



sin 2/? ~ 2- 
A 



-TTItt 



(D.3) 



where 



/= [^M + 7\ — = TTT^ , ■U-4 

J 2^ H^l 2 cos2/? tan2/3-l ' ^ ^ 

and we dropped the A-terms in Eq. (|D.2| ). Notice that, because the term on the left-hand side 
of Eq. (p.l| ) is much smaller then each of the terms on the right-hand side, / ~ A^x^. 

For the purpose of the following estimates we will keep only the largest terms in the variations. 
According to the numbers presented after Eq. (p.4|) , for a small variation of A the largest variation 
on the right-hand side of Eq. ( p.l| ) is 2Ax^AA. We can therefore write 

, , 2Ax2AA 

Mv^)^-A — - — . p.5 

We will justify this approximation a posteriori. Also, in our analysis we will completely neglect 
the dependence of the soft-breaking masses- squared on A and k. This dependence is very weak, 
as seen in the numbers presented after Eq. (|5.4|). 

A small change in A results in a large change in v. Hence, to determine the corresponding 
change in a;, one can use Eq. ( |D.2D and only consider the variation of f ^, which is approximately 
given by Eq. ( p.5|) . We find 

2Ar^ AA 

4fc2xAx ^ \{ksin2(3 - \)^{v'^) ^ A(/tsin2/3 - A)(-4 — - — ) , (D.6) 



9^ 



so that 

Ax _ 2A3(fcsin2/3 - A) AA 
X k'^g'^ A 



(D.7) 



For the point considered in the text (A = 0.11, k = —.045, tan/5 = —2.9) one finds {Ax)/x ~ 
0.2(AA)/A. 

Under a small change Ak, again using Eqs. (p.2| , p.5| ), 

2A^xAx 

AikAkx^ + k'^xAx) ^ X(ksm2(3 - A) (-4 — - — ). (D.8) 

9^ 

Solving for Ax/x. 

Ax Ak ( 2X^iksm2/3-X) Y' 

v^xr^ fcv — j ■ ^ ^ ^ 
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Numerically, {Ax)/x ^ —1.2{Ak)/k. 

Next, we consider the effect of AA on /. The problem comes down to estimating A(cos 2/5)~^, 
which can be done with the aid of Eq. ( |D.3|) : 

1 sin 2/5 . , sin 2/3 / k —m'ir \ , , 



Thus, 



A/ fc(-m2,J sin 2/? AA 



f ' A A2x2 cos3 2/? A ■ ^^'^^^ 

Plugging in the numerical values of the parameters, we find that the right-hand side of Eq. ( p.ll[ ) 
equals — 0.5AA/A. Thus, a 1% change in A results in a 0.5% change in the value of /. Since 
A^x^ changes by 2% in this case, the contribution of / to the variation of v is approximately one 
fourth of that of A^x^, consistent with the numbers given in Subsection 

The above argument can be repeated to find the effect of Ak on /. Notice that sin 2/3 depends 
on the ratio k/X (Eq. ( p.3| )), and hence changing khj +1% has the same effect on / as changing 
A by -1%. 

Finally, we show that the condition for v to remain constant is Ak/k = AA/A. We have 
already argued that sin 2/3, and therefore /, stays unchanged in this case and now show that the 
same is true for A^x^. Under A — A + AA the term A^x^ changes by 2A^x^((AA/A) -|- (Ax/x)) = 
2AV(1+0.2)(AA/A), while under k k+Ak it changes by 2AV(Ax/x) = 2AV(-1.2)(AA;/A;). 
These variations can be made to cancel by imposing Ak/k = AA/A. 
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